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Statistical behavior and scaling properties of iso-height lines in three different saturated two- 
dimensional grown surfaces with controversial universality classes are investigated using ideas from 
Schramm- Loewner evolution (SLE,^). We present some evidence that the iso-height lines in the bal- 
listic deposition (BD), Eden and restricted solid-on-solid (RSOS) models have conformally invariant 
properties all in the same universality class as the self-avoiding random walk (SAW), equivalently 
SLEg/s- This leads to the conclusion that all these discrete growth models fall into the same uni- 
versality class as the Kardar-Parisi-Zhang (KPZ) equation in two dimensions. 

PACS numbers: 05.40.-a, 68.35.Ct, 61.05.-a, 89.75.Da 



Nonequilibrium growth processes exhibit nontrivial 
scaling behavior which are often characterized and clas- 
sified by three exponents, the roughness exponent a, the 
dynamical exponent z and the growth exponent p [1- 
6] . Analytic results for the values of these exponents are 
scarce and one has to depend on numerical analysis. In 
some cases, namely the two-dimensional (2D) Kardar- 
Parisi-Zhang (KPZ) equation, numerical results are not 
definitive and ambiguities remain. On the other hand it 
is by no means clear that this set of exponents is exhaus- 
tive and other characterizing exponents may exist. In 
this paper, we propose a new method for analysis of two 
dimensional rough surfaces based on Schramm-Loewner 
evolution [7]. In this method iso-height lines in the sat- 
urated regime are analyzed as random simple paths in 
which no self-crossing occurs. This leads to an extra 
characteristic for a grown surface, namely /^, the diffu- 
sivity coefficient of SLE. Our approach arrives at a sharp 
conclusion on the universality class of growing surfaces 
in two dimensions. 

The discrete ballistic deposition (BD) [8], Eden, and 
restricted solid-on-solid (RSOS) models [9] (for a review 
and definitions of these models see [2]), are believed to 
be in the same universality class as the KPZ equation [3] 
which describes a non-conserved growth. 
The evidence that two models belong to the same uni- 
versality class can be given in many ways. One of the 
most direct approaches, as pointed out in [10], is to show 
that the two models correspond to the same fixed point 
system. In this direction and based on master-equation 
approach, the KPZ equation has been exactly derived 
for the RSOS model [11], which indicates that these two 
models belong to the same universality class. 
However, the story for the BD model is more contro- 
versial. A continuum equation is derived from the BD 
microscopic rules in [10] which deviates from the KPZ 
equation (the model which is considered in [10] is the 
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next-nearest- neighbor (NNN) BD model, slightly differ- 
ent from the nearest-neighbor (NN) BD model in our 
present paper). Despite this deviation, the symmetry ar- 
guments suggest that the ID BD system is in the same 
universality class of the KPZ equation while for the 2D 
case, the absence of the rotational symmetry in the de- 
rived continuum equation violates the a priori reason for 
them belonging to the same universality class. An ex- 
act lattice Langevin equation for the BD model has been 
derived in [12], whose continuum limit is shown to be 
dominated by the KPZ equation. Although for a ID sub- 
strate the solution of the exact lattice Langevin equation 
yields the KPZ scaling exponents, but for a 2D substrate 
its scaling exponents are again different from those ob- 
tained from simulations [12]. 

Another way to determine the universality class of a 
rough surface is to compute its exponents a, j3 and z. Nu- 
merical results are consistent with the proposition that 
RSOS and KPZ models belong to the same universality 
class, but the situation is more controversial when con- 
sidering the BD and Eden models. 

For the KPZ equation in d = 1, the exact values 
a = 1/2 and /3 = 1/3 are known [3]. The estimated val- 
ues obtained by various numerical works on BD in d = 1 
for roughness exponent and growth exponent range from 
a = 0.42 to 0.506 and /3 = 0.3 to 0.339 [13-17]. Among 
the results, those obtained by Reis [17] are close enough 
to the exact KPZ values. In (i = 2, there is no exact 
computation of the exponents for the KPZ system, nev- 
ertheless, various numerical and theoretical approaches 
have been applied to measure the exponents. The simu- 
lations based on direct numerical integration of the KPZ 
equation in 2D give a = 0.37 to 0.4 [18, 19], and the val- 
ues obtained by various theoretical methods range from 
a = 0.29 to 0.4 [20-23]. Among the theoretical ap- 
proaches, application of the mode-coupling approxima- 
tion for the KPZ equation in 2D [22] yielded a c^ 0.38, 
in good agreement with the values found from simula- 
tions. However, the result a o:^ 0.29 obtained by the 
self-consistent expansion for 2D KPZ equation [20] dis- 
plays a discrepancy with the results of simulations. 



TABLE I: Scaling exponents for BD, Eden, RSOS and KPZ 
[30] models in two dimensions obtained by our simulations. 
With the exception of KPZ, the averages were taken over 10"^ 
independent simulation runs for different square substrates of 
size 50 < L < 700. 



Model 


a 


z 





BD 


0.28(2) 


1.70(5) 


0.15(1) 


Eden 


0.36(2) 


1.65(5) 


0.205(15) 


RSOS 


0.393(10) 


1.58(3) 


0.240(5) 


KPZ 


0.37(1) 


1.61(3) 


0.23(1) 



The growth simulations were undertaken on an 
anisotropic geometry i.e., strips of size L^ x Ly with 
Ly = L, Lj: = 3L and 50 < L < 10^ (for Eden model, 
due to the long CPU time, the simulations were done up 
to size L = 600). Periodic boundary conditions were ap- 
plied in both directions (Strip geometry was chosen for 
congruence with dipolar SLE). We gathered a number of 
5 X 10^ of saturated samples of each size for each of the 
three RSOS, BD and Eden models. The samples were 
obtained from 10^ independent simulation runs. During 
each simulation run, each of the 50 samples was selected 
after 10^ time steps (in units of the number of lattice 
sites) after the saturation time. 



The diversity of the obtained values of various simula- 
tions for 2D BD model ranging from 0.26 to 0.38 for a 
and 0.21 to 0.24 for /3 [14-17, 24, 25], does not indeed pro- 
vide a convincing evidence that it belongs to the same 
universality class of the KPZ model. The same story 
holds for the Eden model with scattered reported results 
for a ranging from 0.20 to 0.39 [2] in (2+1) dimensions. 

We have carried out extensive simulations of the RSOS, 
NN-BD, and Eden models to estimate the values of the 
three exponents a, z and P in (2+1) dimensions. In 
some cases our results take values out of the above men- 
tioned ranges (see table I). Although the two well-known 
scaling relations /3 = a/z and a -\- z = 2 [26] are obeyed, 
within the statistical errors. Nevertheless, due to the sig- 
nificant difference between the exponent values of these 
three models, it is not possible to conclude that they be- 
long to the same universality class. 

A new tool for study of domain walls in critical systems 
is the theory of Schramm-Loewner Evolution (or SLE^) 
[7] (for a review see [27]). The diffusivity constant tz 
determines the critical exponents hence the universality 
class of the system in question. Some authors have ar- 
gued that SLE may be applied to turbulence [28] and 
surface growth phenomena [29-31] as well. Recently, we 
have reported some evidence of conformal invariance in 
the statistical properties of iso-height lines in the satu- 
rated growth models including an experimentally grown 
WOs surface [29], and numerical study of the KPZ equa- 
tion done by direct integration of the discretized KPZ 
equation in (2+1) dimensions [30, 31]. The mere observa- 
tion of scale invariance in a 2D physical system does not 
necessarily imply conformal invariance [32, 33]. However 
arguments pointing to conformal invariance of 2D grow- 
ing surfaces have been attempted [34] . This suggests the 
need for more stringent tests of conformal invariance in 
such systems. 

In this paper, we report the results of extensive sim- 
ulations on three different discrete growth models i.e., 
RSOS, NN-BD and Eden models. We present evidence 
that iso-height lines on saturated surfaces are SLE curves 
of the same diffusivity k, = 8/3 indicating that in this 
spirit and within statistical errors, all these models be- 
long to the KPZ universality class. 




FIG. 1: (Color online) The positive- height clusters shown in 
different colors, and the corresponding spanning iso-height 
lines (solid lines). 



For each height configuration {/li}, a level cut is made 
at the mean height, say h := 0. Then the cluster 
heights were defined as individual sets of connected sites 
with positive height which were identified by using the 
Hoshen-Kopelman algorithm [35]. All spanning clusters 
along ^—direction were marked. For each spanning clus- 
ter, a walker algorithm was applied to determine each 
of its perimeters which connects the lower boundary to 
the upper one i.e., the spanning iso-height lines (Fig. 1). 
The iso-height lines were identified uniquely by using the 
tie-breaking rule on the square lattice, described in [36]. 
Thereby, an ensemble of contour lines of fixed linear size 
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FIG. 2: (Color online) The average length / of a spanning iso- 
height line of the saturated growth models versus the system 
size L, I ^ L^f . The solid line shows the expected result for 
SAW. The error bars are almost the same size as the symbols. 



L was obtained for each surface ensemble of different size. 

We base our arguments on four different tests acertain- 
ing that the iso- height hues are SLE: the fractal dimen- 
sion, winding angle statistics, left passage probability and 
direct SLE test. 

i) Fractal dimension. For conformally invariant curves 
the fractal dimension is related to the diffusivity n by 
the relation dj = 1 + n/S. The fractal dimension and the 
conjectured value of diffusivity for SAW are known to be 
df = 4/3 = 1.33 and /^ = 8/3 = 2.66. 

In Fig. 2 we show our computed results for the frac- 
tal dimension of the iso- height lines for different models. 
The best fit to our data in the whole range of the lattice 
sizes yields df = 1.345(7), 1.330(5) and 1.335(4) for BD, 
Eden and RSOS models, respectively. In all of our mea- 
surements for BD model, our analysis on slightly larger 
system sizes 150 < L < 10^ corresponds to the same re- 
sults as the Eden and RSOS models. For example we 
find df = 1.337(5) for the BD model in this size range. 
As shown in Fig. 2, our results are well compatible with 
the fractal dimension of SAW and SLEg/s. 

ii) Winding angle statistics. The winding angle be- 
tween two end points of a finite SAW in two dimensions, 
is studied in [37] using Coulomb gas methods. They 
found that the winding angle is Gaussian distributed with 
a variance of ^ (8/^)lnL, where L is the distance be- 
tween the end points and g is Coulomb gas coupling pa- 
rameter which is related to k, hy g = A/k,. They have 
also shown that the winding angle at a single end point 
relative to the global average direction of the curve is a 
Gaussian with variance of (4/^) InL. 
It is shown in [38] that the variance in the winding {0'^) 
at typical points along the curve is 1/4 as large as the 
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FIG. 3: (Color online) Variance of the winding angle for the 
spanning iso-height lines of the saturated growth models. The 
solid lines are set according to the Eq. (1) for comparison, 
with appropriate obtained values of a for each model and k, = 
8/3. In the inset, the variance in semilogarithmic coordinates. 
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FIG. 4: (Color online) Left passage probability computed for 
the spanning iso-height lines of different models with p — 
O.IL. The solid line shows the prediction of SLE for n = 8/3. 



variance in the winding at the end points, i.e.. 



{0^ 



(K/4)lnL. 



(1) 



Using the same definition as in [38], we measured {6'^) for 
different models, results are shown in Fig. 3. Data points 
compare well with Eq. (1) (solid lines), using tz = 8/3 
and a suitable value of the parameter a obtained from the 
best fit to data for each model. The direct measurement 
of K: also obtained from the best fits to the data shown 
in the inset of the Fig. 3. We find almost the same value 
tv/A = 0.700(20) for BD (again for larger sizes), Eden and 
RSOS models, within the statistical error. 

iii) Left passage probability. The probability P^i^) that 
an SLE^ curve, in the upper half-plane, passes to the left 
of a given point at polar coordinates (p^^p)^ is computed 
by Schramm [39] 



1 r(i) 



(2) 

where 2^1 is the hypergeometric function. 
As another check, we measure this quantity (which 
should also hold for a dipolar SLE near the starting point 
i.e., p <C L [40]) for the contour lines of different growth 
models. 

As shown in Fig. 4, our results for the three models are 
again in good agreement with the prediction for SLEg/a. 
iv) Direct SLE test. Using a discrete Loewner evolution 
and successive appropriate conformal maps, we extracted 
the Loewner deriving function {^t} of each iso-height line 
represented by the sequences of points {zq, ^1, ..., ^at} in 
the complex half-plane, with zq = (0,0). We use the 
function gt{z) = ^{z — ^t) + 4t + ^t, with t = \^zf and 
£,t = 3?^i , to map all of the points except the first one to 
a shortened renumbered sequence. After each recursive 
map the first point in the sequence is swallowed and a 
sequence of {^^} can be obtained for each iso-height line. 
We have also checked the map appropriate for dipolar 
SLE [40], and found no significant differences. 
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FIG. 5: (Color online) The diffusivity k obtained for each 
contour ensemble versus their linear size L. Inset: statistics 
of the driving function {^t } obtained for the contour ensemble 
of the RSOS model with L = 10^ 



We find that the statistics of the deriving function for 
each curve ensemble of fixed linear size L, converges to 
a Gaussian process with variance (^|) = f<it. Finite size 
effects may be reduced by looking at shorter segments of 
the curve e.g., when 10% of the total average length of 
the curves is mapped. An example is shown in the inset 
of Fig. 5, obtained for the contour ensemble of the RSOS 
model with L = 10^. We observe that hi shows a slight 
dependence on system size, but reducing with L. We find 
that the value of the diffusivity k, for larger system sizes 
approaches the expected value for SAW i.e., k = 8/3, for 



all three models. As shown in Fig. 5, for RSOS model 
this convergence begins in rather smaller sizes but for 
Eden and BD models larger system sizes are needed. 

The more accurate results obtained here for the BD 
model is slightly different from that reported in [29] , this 
is due to the considerable difference in the number of 
averaging samples and reduced finite size effects. 

Summing up, although the numerical values of three 
exponents a, /3 and z for BD, Eden, RSOS and KPZ 
models are scattered, the numerical value found for Hi is 
sufficiently sharp to suggest that these models all belong 
to the same universality class. What remains is the inter- 
dependence of these exponents and tz. The existence of 
two scaling relations /3 = a/z and a -\- z = 2^ guarantees 
that there is only one independent exponent e.g., the 
roughness exponent a. On the other hand, this paper 
introduces a new characteristic value for a rough surface: 
hi. How inter-dependent are these two? There is a pow- 
erful scaling argument given in [41] which connects df of 
a contour line to a of the same surface, df = 2 — xi — a/2, 
where xi is the loop correlation exponent. Although the 
exact value of xi = 1/2 is known only for the limiting 
cases of a = and 1, but it is conjectured that its value 
is super universal and is independent of a for Gaussian 
surfaces. This leaves the case of 2D KPZ in ambiguity, 
since it does not follow a Gaussian distribution. A simple 
minded value xi = 1/2 leads to df = 4/3 giving a = 1/3 
which although elegant [21] but is in conflict with nu- 
merical results. Our future efforts will concentrate on 
revealing the nature of this relationship. 
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